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Abstract. We investigate the AC conductivity of binary random impedance networks, 
with emphasis on its dependence on the ratio h = ai/ao, with (Tq and ai being the 
complex conductances of both phases, occurring with respective probabilities p and 1 — p. 
We propose an algorithm to determine the rational /i-dependence of the conductance of a 
finite network, in terms of its poles and of the associated residues. The poles, which lie on 
the negative real /i-axis, are called resonances, since they show up as narrow resonances 
in the AC conductance of the RL — C model of a metal-dielectric composite with a high 
quality factor Q. This approach is an extension of a previous work devoted to the dielectric 
resonances of isolated finite clusters. A numerical implementation of the algorithm, on the 
example of the square lattice, allows a detailed investigation of the resonant dielectric 
response of the binary model, including the p-dependence of the density of resonances 
and the associated spectral function, the Lifshitz behaviour of these quantities near the 
endpoints of the spectrum of resonances, the distribution of spacings between neighbouring 
resonances, and the Q-dependence of the fraction of visible resonances in the RL — C model. 
The distribution of the local electric fields at resonance is found to be multifractal. This 
result is put in perspective with the giant surface-enhanced Raman scattering observed 
e.g. in semicontinuous metal films. 
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1 INTRODUCTION 



Random networks of complex impedances are currently used to model electrical and 
optical properties of disordered inhomogeneous media. The most common situation is that 
of a binary composite medium, modeled by attributing a random conductance to each bond 
(x, y) of a lattice, according to the binary law 

{(To with probability p, (II) 
ui with probability q = 1 — p, y • ■> 

in correspondence with the bond percolation problem (see ref. [1] for a review). 

The conductances (inverse impedances, or admittances) (Tq and ai of both phases 
take arbitrary frequency-dependent complex values. Hereafter we follow the notations of 
ref. [2] . The dimensionless complex ratio 

h=^ (1.2) 

and the concentration p are the essential parameters of the model. As far as static (DC) 
properties are concerned, the limiting case h = embraces the conductor-insulator mixture 
(o"! = 0) and the superconductor-conductor mixture (ctq = oo). In both situations the 
conductivity exhibits power-law behaviour for p close to the critical concentration p^, 
corresponding to the percolation threshold. 

More generally, the conductivity of the binary model obeys a scaling law in the critical 
regions defined by |/i| ^ 1 and |i? — Pc| ^ 1, or equivalently ^ 1 and \q — Pc\ <S 1 [!]• 
Frequency-dependent (AC) properties of metal-dielectric composites, such as cermets or 
thin films, are often investigated by means of either the R—C and the RL—C models [3]. In 
both cases the low-frequency regime corresponds to \h\ ^ 1, and thus to critical behaviour 
when the metallic concentration q is near the percolation threshold. 

The purpose of the present work is to shed some new light on the resonant behaviour 
of the binary model. The emphasis will be put on the analytic structure of the conductivity 
in the ratio h, or in the equivalent complex variable 

X = -l^ = ^^. (1.3) 
1 — n (To — cTi 

For any value of the concentration p, the conductivity has singularities for h real and 
negative, i.e., in the range < A < 1. A quantitative investigation of these singularities, 
and of related quantities, is the main goal of the present paper. Our aim is twofold. First, 
the analytic structure of the conductivity of the binary model is a classical subject, since 
the developments of the Bergman-Milton theory [4, 5]. The key ingredient of this formalism 
is the spectral function H{p,x), which has only been the subject of a limited number of 
investigations so far [6, 7]. The present approach provides a direct accurate numerical 
evaluation of the spectral function of binary random networks. Second, the singularities 
of the conductivity have a physical interpretation in terms of dielectric resonances in the 
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RL — C model, and of the relaxation times in the transient response in the R — C model. 

The regime of dielectric resonances has been argued to provide a natural explanation for 
the anomalous fluctuations of the local electric field [8], which are responsible for giant 
surface-enhanced Raman scattering observed e.g. in semicontinuous metal films [9]. 

The dielectric resonances of isolated clusters have been investigated in ref. [2] . The sit- 
uation considered there was a finite set of (metallic) bonds with conductance ai , embedded 
in an infinite (dielectric) host lattice, whose bonds have a conductance ctq- It was shown 
there that the conductance of such a system is entirely characterised by a finite number of 
resonances. The positions of the resonances, and the associated cross-sections 7^, are 
expressed in terms of the eigenvalues and eigenvectors of a finite matrix M, dictated by 
the geometry of the clusters. The present paper is an extension of the method of ref. [2] 
to an arbitrary binary network. We shall make use of an eflBcient algorithm, which allows 
for an exact determination of all the resonances of a finite sample. The setup of this paper 
is as follows. In section 2, we gather definitions and results on various features of the con- 
ductivity of the binary model. Section 3 is devoted to the presentation of the algorithm. 
Section 4 contains a variety of numerical results, concerning especially the spectral func- 
tion of the Bergman-Milton theory and the spectral density of resonances, their Lifshitz 
behaviour near the endpoints of the spectrum (A = and A = 1), the number of visible 
resonances of a finite sample and its dependence on the quality factor Q of the RL — C 
model, and the distribution of the local electric fields at resonance, which turns out to be 
multifractal. 

2 THEORETICAL BACKGROUND 

2.1 Conductance and impedance of a finite network 

Consider the binary model on a finite network of size M x A^, as shown in Figure 1. 
This network contains = (M - 1)N sites (nodes) and hb = MN (M - l)(iV - 1) 
bonds (links), among which MN are horizontal (perpendicular to the electrodes), and 
(M — 1)(A'" — 1) are vertical (parallel to the electrodes). The bonds with a conductance uq 
are called the P-bonds, which form the V-set. There are n-p of them, among which are 
horizontal and rij, are vertical. Similarly, the bonds with a conductance ai are called the 
Q-bonds, which form the Q-set. There are ng of them, among which ng are horizontal 
and rig are vertical. 

Let Y be the conductance (admittance) of the network, measured between the plane 
electrodes (bus bars) shown in Figure 1, and let Z = 1/y be its impedance. These 
quantities are rational functions of the dimensionless complex variables h or A. Let us 
anticipate that it is more convenient to use the variable A. We have 




(2.1a) 



The prefactors of these expressions are the conductance and the impedance of the uniform 
network whose bonds have a conductance ctq, since A = oo corresponds to a"i = ctq (no 
disorder). The involved in the product expressions are the poles of the conductance 
and the zeros of the impedance, while the Aq are the zeros of the conductance and the 
poles of the impedance. The number nn of poles and zeros depends on the configuration 
of the random bonds. It is bounded by the number of sites of the network: < nu < ns- 
Finally, the partial- fraction expansions involve residues aa and Pa- 

It follows from considerations on the dissipated power that Y/ao is a Stieltjes function, 
namely its imaginary part has the sign of Im/i, or equivalently of ImA [4, 5, 1]. This 
property implies that the poles and zeros of Y and of Z alternate, according to 

< Ai < Ai < • • • < A„^ < A^^ < 1. (2.2) 

We have Ai = if, and only if, the Q-set is conducting, i.e., it contains at least one 
connected path between both electrodes. The conductance of this Q-set then reads 

N 

Yq = (2.3) 

Similarly, we have A^^j = 1 if, and only if, the V-set is not conducting. 

The Stieltjes property also implies that the residues and (3a are positive. These two 
sets of residues are different from each other in general. They obey, however, a remarkable 
sum rule, which can be proved by expanding the conductance Y around A = oo, to first 
order in 1/A. Indeed, 1/A = (ctq — cri)/ao is the dimensionless contrast between the 
conductances of both phases. As a consequence, for a given realisation of the random 
network, up to first order in 1/A included, the conductance reads Y = {N/M)aH, where 

^ MN "V MNX) ^ ' 

is the average value of the conductances of the horizontal bonds of the network. By 
inserting this estimate into the second expressions of eqs. (2.1a, lb), we obtain the relation 

nR riR H 



In the case of isolated finite clusters [2] , the poles and zeros of the conductance form 
very tight doublets. We have indeed 

aa-^Pa-^Xa-Xa-^ " ' (^.6) 

The sum rule (2.5) then becomes 

Y,laXa{l-Xa)^nH, (2.7) 



a=l 
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where uh is the number of horizontal bonds in the clusters (perpendicular to the elec- 
trodes) . This identity was not noticed in ref . [2] . 



2.2 R-C model 

The R — C model, already mentioned in the Introduction, is defined as follows. The 
Q-set is a metallic phase, whose bonds consist of a pure resistance R, while the V-set is a di- 
electric phase, whose bonds consist of a perfect capacitance C. The complex conductances 
at frequency / = uj/{2'k) thus read 

(jQ = iCu), (71 = 4- (2-8) 

Jri 

Along the lines of refs. [1, 2], we introduce the microscopic relaxation time 

T = RC, (2.9) 



so that 



h=^, X=-^. (2.10) 
iujT iujT — 1 



We thus have /j. — > oo and A — > at low frequency. 

The poles of the conductance Y show up as relaxation times in the transient response 
of the model [1]. Consider indeed the R — C model on a finite network, submitted to a 

voltage V{t) = Vo9{t) applied between the electrodes, with 9{t) being the Heaviside step- 
function. The intensity I{t) across the network in the transient regime can be evaluated 
by means of the Fourier transformation. We thus obtain 

7_oo 27r za; + 

The second expression of eq. (2.1a) for the conductance Y, with A given in eq. (2.10), 
yields the result 

where the relaxation times read 

Ta = ^^T. (2.13) 

If the metallic Q-phase is conducting, we have Ai = 0. The term a = 1 in eq. (2.12) 
yields the DC current through the network, Iq = NVoai/{MR), in agreement with the 
result (2.3). 

2.3 RL-C model 

In the RL — C model, already mentioned in the Introduction, the metallic bonds of 
the Q-set now consist of an inductance L in series with a weak resistance i?, while the 
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dielectric bonds of the T'-set still consist of a perfect capacitance C. The conductances at 
frequency f = u)/{2'k) now read 

(To = iCcj, ai = / ■ (2.14) 

it + zLo; 



Along the lines of refs. [1-3] , we introduce the microscopic resonance frequency 

1 

7W' 



the reduced frequency 
and the quality factor 



u> 



y=—, (2.16) 



which is a dimcnsionlcss measure of the dissipation rate. 

In the following, we shall mostly consider the case of a weak dissipation, corresponding 
to a large quality factor {Q 3> 1). We have then 

7 _ 1 1 i y'^ — iy/Q iy /o ■\fi\ 

^ -y^ + iy/Q ^~¥~¥Q' ^ 1 + - iy/Q ^ i+7 ~ {i + yYQ' ^ ' ^ 

We notice that the low-frequency regime {y — > 0) again corresponds to /j. — > oo and A — > 0. 

Since the variables h and A have a small (negative) imaginary part, proportional to 
1/Q, the poles of the conductance Y and of the impedance Z of the network show up in 
the frequency dependence of these quantities as narrow resonances, respectively located at 
<^a = i^oya and uja = i^oya, with 



We have indeed 



Y 



NCiVQ ^ a, l/{2Q)-t{y-ya) 
2M ^(1-A„)2 {y-y,r + l/{4Q^y ^ ' ^ 



M ^ Pa l/{2Q)-i{y-y^, 

" - 2NCu^, ^ A„(i - K) (y - ya? + i/(4Q2) • ^ • ^ 

Eq. (2.20a) shows that the real part ReF(a;) of the conductance exhibits nn narrow 
resonances, at Ua = <^o2/o- The resonance peaks have a Lorentzian shape, with a common 
absolute width Aa; = ujq/{2Q). The maxima at resonance read 

(ReF) (2.21) 
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while the area under each resonance peak is 

A similar pattern of resonances can be observed on the real part Re Z(lj) of the impedan- 
ce (2.20b), with maxima at resonance and areas under resonance peaks respectively given 

by 

. X ML Pa J TtM Pa 

2.4 Spectral function and density of resonances 

The conductance Y of an infinitely large network drawn on any regular lattice becomes 
a self- averaging quantity. In other words, the conductivity 

MY 

S(j9,A)= lim — — (2.24) 

is an intrinsic characteristic of the binary model, which depends on the lattice under 
consideration, on the concentration p, and on the complex variable h or A. The conductivity 
of the binary model has been the subject of many investigations [1]. We gather below some 
definitions and properties which will be useful in the sequel. 

The second expression of eq. (2.1a) yields the Bergman-Milton integral representa- 
tion [4, 5] 

.„(,_^'£|^), (...5) 

where the positive function 

H{p,x)^ lim \2aa6{x-Xa) (2.26) 

M,iV— »oo 

a=l 

is called the spectral function of the binary model. In more precise terms, it is the density 
of a positive measure, supported by the interval < a; < 1. Indeed, as a consequence of 
eq. (2.3), the spectral function has a singular component at a; = 0, of the form 

H,^{p,x) = A{l-p)5{x), (2.27) 

with the notation to be introduced in eq. (2.38), whenever the Q-set is conducting, i.e., 
for 1 — p = q > pc- 

The spectral function entirely determines the conductivity S(p, A), by means of the 
integral representation (2.25). Conversely, it is given by the inverse formula 

1 ^ X + iO) , , 

H{p,x) = -Im^ L. (2.28) 

TT (Jo 



In other words, the conductivity is analytic in the complex A-plane cut along the interval 
< A < 1, and its discontinuity along this cut reads 

Disc S(j9, A) = 2i7raoH{p, A). (2.29) 



The spectral function directly yields the transient intensity I{t) of the R — C model 
on a very large network, namely 



To close up, we introduce the spectral density of resonances 
and the total density of resonances 

^^^^^ = M,feoo L ^^-^^^ 

representing the mean number of resonances per site. This quantity has a non-trivial 
dependence on the concentration p. It will be expressed in eq. (3.14) in terms of geometrical 
quantities of the bond percolation problem. 

2.5 Homogeneity and duality 

First, the conductivity of the binary model is invariant under the simultaneous in- 
terchange p q = 1 — p, (To (Ti. Indeed, let Y = aoF{h) = (To/(A) be the con- 
ductance of the network G, shown in Figure 1. Then the network G', obtained from G 
by interchanging all the bond conductances according to ctq ^ cri, has a conductance 
Y' = aiF{l/h) = (7i/(l — A). For a large enough network, if G is typical of the concentra- 
tion p, then G' is typical of the concentration 1 —p. The conductivity, the spectral function, 
and the density of resonances therefore obey the following identities, on any regular lattice 



E(p,/i) = /iE(l -p, l//i), (2.33a) 

AS(p,A) = (A-l)S(l-p, 1-A), (2.33b) 

xH{p, x) = (1 - x)H{l - p, 1 - x), (2.33c) 

p{p,x) — p{l — p,l — x). (2.33d) 



More interestingly, the square lattice is self-dual, i.e., invariant under the geometric 
transformation called duality. This concept has been introduced in physics by Kramers and 
Wannier [10], while its consequences on random resistor networks have been explored in a 
systematic way by Straley [11]. The dual G of a planar network G has bonds which cross 
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those of G, the conductances of any pair of crossing bonds being inverse to each other. The 
duahty property imphes that the conductance of the whole network G reads Y = 1/Y = Z . 
With the same notations as above, we have Y = {l/ao)F{l/h) = (1/(To)/(A), so that 
F{h)F{l/h) = /(A)/(l - A) = 1. Therefore 

S(p,A)S(p,l-A) = (72, (2.34) 

and 

p{p,x) = p{p,l-x). (2.35) 

Because of its non-hnearity, the duahty identity (2.34) for the conductivity does not yield 
any identity involving the spectral function H{p, x) only. 

The above identities can be combined with the homogeneity properties (2.33). We 
thus obtain that the percolation threshold is — 1/2, and that the conductivity and the 
spectral function right at this point read 



E(l/2, A) = ctoa/^, H{1/2,x) = - (2.36) 

V A TT \ X 

Furthermore, the spectral and total densities of resonances obey the relations 

p{p,x) = p(p,l- x) = p{l-p,x) = p{l-p,l- x), pr{p) = pr{1-p). (2.37) 

2.6 Critical behaviour and scaling 

As mentioned in the Introduction, the conductivity exhibits scaling behaviour around 
the percolation threshold p — pc- The conductivity of the conductor- insulator mixture 
vanishes for p <Pc, while it reads 

S((7i = 0) = aoA{p) {p > Pc). (2.38) 

Similarly, the conductivity of the superconductor-conductor mixture is infinite for p > Pc, 
while it reads 

S((7o = oo) = aiB{p) {p < Pc). (2.39) 
Both amplitudes A(p) and B(p) have a power-law behaviour for p near pc'. 

A{p) ^ a{p - pc)^ (p^pt), /2 40) 

B{p) f^abipc-p)'" ip^Pc)- 

The conductivity of the binary mixture obeys a scaling law in the critical region 
defined by <^ 1 and \p — Pc\ <S 1, of the form [1] 

S(p, h) ^ao\p- Pel' *± (h\p - Pel"'"') , (2.41) 
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where the ^± are scahng functions of one complex variable, with ± referring to the sign 
of p — Pc- The homogeneity relation (2.33a) allows to describe the vicinity of the other 
critical point, \h\ 3> 1 and |g — Pc\ <^ 1- 

The scaling formula (2.41) reproduces the power laws (2.40) for small values of the 
argument of the scaling functions, as we have a — $_(_(0), while $_(a;) ~ 6x as x — 0. On 
the other hand, both scaling functions have the common power-law behaviour ^±{x) ~ 
Kx'^, with u = t/{s + t), as |a;| — > oo and |Arga;| < tt, hence 

S(pe, h) Kao h"" fti K{a'oa{) (2.42) 

for |cri| ^ |o"o|, right at the percolation threshold. More generally, the power-law be- 
haviour (2.42) holds for hi, <C \h\ <C 1, where the crossover scale /i* reads 

Kr.\p-Pc\'+\ (2.43) 

while the scaling laws (2.38-40) hold in the opposite regime {\h\ ^ /i*). 

As a consequence of eq. (2.28), the spectral function also obeys a scaling law of the 
form (2.41), namely 

H{p,x)^\p-Pc\' F±(^{l-x)\p-pc\-'-'^, (2.44) 

for 1 — X <^ 1 and |p — Pc| -C 1, and a similar law around the other critical point, i.e., for 
X <^ 1 and \q — Pc\ <^ 1- We shall come back to the scaling law (2.44) in section 2.9. 

Consider now the binary model on a large but finite sample, of size M x N. In the 
critical region, its mean conductance obeys the finite-size scaling law 

Y « aoN-^/" ^ (^{p - Pc)N^'\ {p - Pc)h-^'^'^*\M/N^ , (2.45) 

where \1/ is a three-variable scaling function. As a consequence, right at the percolation 
threshold, the critical region extends over a range 

5h ~ Ar-(«+*)/^. (2.46) 

On the square lattice, the duality symmetry implies Pc = 1/2, and A{p)B{l —p) = l, 
hence s — t and u = 1/2, in agreement with eq. (2.36). The common numerical value of 
these exponents is s/v = t/u = 0.9745 ± 0.0015 [12], with the exponent of the correlation 
length being exactly u = 4/3, hence s = t = 1.300. Hence, for p = Pc = 1/2, the same 
critical singularity simultaneously affects both endpoints of the spectrum, A = 1 and A = 0, 
corresponding respectively to h = and h — oo, again in agreement with eq. (2.36). 

2.7 Sum rules 

The representation (2.25) of the conductivity can be expanded as the following series 
in inverse powers of A 

E(p,A) = .„fl-f:^y (2.47) 

\ k=0 J 
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where the coefficients 



Hk{p)= x^H{p,x)dx= Jrm ^c^aA^ (2.48) 

/n M,iV — »oo 

are the moments of the spectral function H{p, x). 

As we have aheady noticed in section 2.1, the expansion variable 1/A is the dimen- 
sionless contrast between the conductances of both phases. As a consequence, the expan- 
sion (2.47) can be viewed as a special case of the weak-disorder expansion of the conductiv- 
ity of a random network with an arbitrary narrow distribution of bond conductances [13]. 
The general results to sixth order derived there can be transcribed in the present case of 
binary disorder on the square lattice. We thus obtain the following expressions for the first 
six moments, 

Ho{p)^q. H,{p)^^, if2(p) = f, 

Hs{p) = ^(1 +pq), H^p) = ^(l + 3pq+{J- l)pq{p - q)), (2.49) 
H^ip) = (1 + ("^ + 5)P? + 4( J - l)pq{p - g) - 2(2 J - 3)^2^2) ^ 

with q=l-p. The number J = = J2 = 1.092958179, with the notations of ref. [13], is 
the only non-trivial quantity occurring in the sixth-order expansion. We have thus derived 
explicit sum rules for the spectral function, which agree with the expressions given in 
ref. [7]. The number J, denoted there as 1 — a^, was estimated there from numerical data 
to be J = 0.9 ±0.5. 

2.8 Effective-medium approximation 

The effective- medium approximation (EMA), introduced by Bruggeman in 1935 [14], 
is a self-consistent approximate scheme to evaluate the conductivity of random impedance 
networks [15, 16, 13], which is still being very widely used [17]. 

In the present case of the binary model on the square lattice, the EMA prediction for 
the conductivity is given by 

yEMA _ yEMA _ 



hence 



.EMA 



(2.51) 



ao ((p - 1/2)(1 -h) + ^(p - 1/2)^(1 -hr + h) 
= ^(p-l/2 + ^{p - 1/2)2 + x{\ - 1)) . 
This EMA formula for the conductivity is analytic in the A-plane cut along the interval 

Amin < A < Amax, with 

Amin = 1/2 - VKI-P), Vax = 1 - A^nin = 1/2 + y/p{l-p). (2.52) 
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The prediction for the associated spectral function [cf. eq. (2.28)] reads 



H 



■EMA 



'max 



X){X - Amin) 



(2.53) 



TTX 



The EMA formula gives a very accurate approximation to the conductivity of the 
binary model in generic circumstances. For instance, the 1/A-expansion of the EMA for- 
mula (2.51) gives expressions Hf^^{p) for the moments which only differ from the true 
results (2.49), starting with H4{p), by replacing J by J^ma ^ ^ ^^gj rpj^^ EMA scheme 
also respects the duality symmetry (2.34). The predictions (2.51), (2.53) therefore agree 
with the exact results (2.36) for the conductivity and the spectral function at p — pc — 1/2. 
The EMA also correctly predicts the equality s = t, but not the common value of these 
exponents {s = t = 1 instead of 1.300). 

For a generic value of the concentration [p ^ pc), the endpoints (2.52) are such that 
< Aniin < Amax < 1- The support of the EMA prediction (2.53) for the spectral function 
thus docs not extend over the whole interval < A < 1. In the vicinity of the percolation 
threshold {p ^ Pc — 1/2), we have Amin = 1 — Amax ~ (p — Pc)'^, in agreement with the 
estimate (2.43) of the crossover scale h^,, with s + t = 2. 

2.9 Lifshitz tails 

It has been argued [6] that the spectral function H{p, x) of the true conductivity of 
the binary model extends over the whole allowed spectrum < A < 1, for any value of the 
concentration p, at variance with the EMA formula (2.53). 

This argument has then been put in perspective [1] with Lifshitz singularities [18-20]. 
These singularities are caused by the presence of very large, and thus very improbable, 
ordered regions in a randomly disordered system. The original example considered by 
Lifshitz [18] is that of the phonon spectrum of a binary harmonic alloy, consisting of 
light atoms, with mass m and concentration p, and heavy atoms, with mass M > m and 
concentration q = 1—p. Lifshitz has argued that the vicinity of the upper edge Wmax of the 
phonon spectrum of the alloy is dominated by large ordered regions, almost spherical in 
shape, consisting only of light atoms. He thus showed that Wmax coincides with the upper 
edge of the pure lattice consisting only of light atoms, and that the density of states of the 
alloy vanishes exponentially fast near a;niax) as 



where c is a lattice-dependent constant, which can be evaluated exactly. Along this line 
of thought, it has been argued in ref. [1] that the spectral function of the binary model 
has an exponentially small Lifshitz tail, extending all the way to the endpoints A = and 
A = 1 of the spectrum of resonances, of the form 




(2.54) 




{x 




(2.55) 
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and similarly for x — >^ 1. This expression was rather conjectural, as the determination of 
the relevant ordered regions was left as an open question, so that the prefactor C{p) was 
not predicted. 

Hesselbo [21] then argued that the relevant ordered regions are hairpin configurations, 
as shown in Figure 2a. Let be the transversal conductance of the hairpin consisting of n 
cells, measured between the point electrodes shown in Figure 2a, considered as an isolated 
network (not embedded in the square lattice). This quantity can be evaluated from the 
recursion relation ^ 

Yn = (Ti+(— + , (2.56) 



with Yq = (7q. Consider first a value of h not on the negative real axis, and set 

/i= — = 2sinhV, (2.57) 

with Re// > and |Im/i| < tt. The Mobius map involved in the recursion (2.56) has two 
fixed points, 

Y^ = ^{e^'^-1), (2.58) 

with the stable fixed point 1+ being the transversal conductance of the infinitely long 
hairpin (ladder). The recursion (2.56) can be solved explicitly, along the lines of ref. [1], 
by means of the variable t„ = (F^ — Y_)/{Yn — i + )- We thus obtain 

3cosh((2n + l)u) — cosh((2n — l)jLi) 
Yn = cTosinh/x }\ . 2.59 

3smh(2n//) — smh(2(n — 1)/*) 

For a strong dielectric contrast (/i ^ 0), the correlation length of currents along the 
hairpin diverges according to ^ = 1/(2/^) ~ (2/i)^^/^. In this regime, the conductance of 
long hairpins scales as 

Y^p^ — 2;cotanh2;, with z = ^ = 2n//. (2.60) 
2n 4 

This scaling form of the conductance exhibits an infinite array of alternating zeros, lying 
at 

,e.. (2.61) 



and poles, lying at 



Za = aiTT, i.e., 1- Xa^ ~2n^' {'^■^'^) 



with a > 1. The first zero Ai, corresponding to a pole in the dual configuration, yields, 
according to Hesselbo, the Lifshitz behaviour of the spectral function, at least for a small 
enough concentration p. 

It turns out that the formula (2.59) also gives the conductance of worm-like networks, 
such as the configuration shown in Figure 2b, where the square cells are put together in 
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any random fashion, respecting the hnear structure and the constraint of self- avoidance. 
The number of such worm-hke configurations with n cells is of order exp(?i5'), with S 
being the associated configurational entropy. On the other hand, a hairpin with n cells 
occurs with a probability of order p'^"', at least for p small enough. Altogether, the first 
zero of eq. (2.61) is expected to show up with a probability weight of order (p^e'^)". By 
eliminating the number n between the above estimates, we obtain the following analytical 
form for the Lifshitz tail of the density of resonances and of the spectral function 

pip,x)^H{p,x)r^exp(^-^^ (x^O), (2.63) 

and a similar formula for x ^ 1. This result, with an inverse-square-root behaviour in 
the spectral variable x, is characteristic of Lifshitz tails in one-dimensional systems [19]. 
This is due to the fact that the relevant structures are linear objects. Furthermore, the 
identities (2.33), (2.37) imply C{p) — C(l — p). The above argument also leads to a 
prediction for the small-p behaviour of the amplitude C{p), namely 

C{p) ^ ^ (^\lnp\ - I) {p « 1). (2.64) 

These predictions will be compared to numerical data in section 4. 

The Lifshitz tail of the spectral function manifests itself in the long-time tail of the 
transient intensity response of the R — C model. Indeed, in the regime of long times 
(t >> r), eq. (2.30) is dominated by the vicinity of a; = 0, hence 
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In the critical region, the prediction (2.63) for the Lifshitz tail is compatible with the 
scaling law (2.44) for the spectral function, provided the amplitude C{p) vanishes near the 
percolation threshold, according to 

C{p) r^\p-Pc\^'+'^/''. (2.66) 

The Lifshitz behaviour (2.63) is expected to hold only deep in the tails, for x (or 1 — x) 
much smaller than the crossover scale /i^, defined in eq. (2.43). 

3 ALGORITHM 

Wc now turn to the presentation of our algorithm for evaluating the rational h- 
dependence of the conductance of a finite binary network, such as that shown in Figure 1. 
This approach is an extension of the method of ref. [2]. It turns out that a very similar 
formalism was proposed by Straley [6] some twenty years ago, but this work has appar- 
ently not been noticed since then. The conductance will be determined in the second form 
of eq. (2.1a), namely we shall calculate first the poles Xa of the conductance, giving the 
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positions of the resonances, and then the associated residues eta, giving the strengths of 
the resonances. 



3.1 Generalities 

Along the hnes of ref. [2], our starting point is the Kirchhoff equations for V^, the 
electric potential at site x: 

J]ax,y(Fx-Fy)=0. (3.1) 

y(x) 

There is one such equation per site x inside the network. The notation y(x) means that 
y is a neighbour of x, i.e., there exists a bond (x, y), and the sum possibly includes sites 
y belonging to either electrode. Eqs. (3.1) have to be complemented by the boundary 
conditions Vy = for the sites y on the left electrode, and Vy = Vq for the sites y on the 
right electrode. 

We define the topological Laplace operator A on the network as 

(AV),= ^(yy-V,), (3.2) 
y(x) 

again with the convention that the sum possibly includes sites y belonging to either elec- 
trode, in which case the Dirichlet boundary condition V"y = is assumed. The operator A 
can be written as the sum A = A-p + Aq of its components A-p on the V-set and Aq on 
the Q-set, respectively defined as 

iA-pV)^= J2 (Vy-V^l (AsF)x= Yl (^y-^x), (3.3) 
yGPCx) yeQ(x) 

where y e P(x) (respectively, y e Q(x)) means that (x, y) is a V-hond (respectively, a 
Q-bond). We also introduce the quantities 



Ax^^ = 1 iff X is a neighbour of the left electrode, 
= 1 iff X is a neighbour of the right electrode, 
= 1 iff X is connected to the left electrode by a Q-bond, 



(3.4) 



B^^ = 1 iff X is connected to the right electrode by a Q-bond. 
With these notations, the Kirchhoff equations (3.1) read 

ao{ApV)^ + ai{AQV)^ + Fo(^o(4''^ - B^J'^) + = 0, (3.5) 

or equivalently, in vector and matrix notation, 

(Aq-AA)V = Vb(AA(^) -B(^)). (3.6) 
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This reduced form only involves the complex variable A defined in eq. (1.3). The conduc- 
tance of the network is given by 

y = ^, (3.7) 

where the total current / flowing into the network from the left electrode reads 

I = E(^o(AW _ + (3.8) 



or equivalently, in vector notation, 

/=^(AAW-B(^))-V. (3.9) 



3.2 Poles of the conductance 

In analogy with ref. [2] , the poles of the conductance are the non-trivial values Xa of A 
for which the homogeneous Kirchhoff equations (3.1) with Vq = have a non-zero solution, 
namely 

(Aq - A«A)V = 0. (3.10) 

This is a well-posed generalised eigenvalue problem, since A and Ag are two real symmetric 
matrices, of size ns x ns, and (—A) is a positive definite matrix. 

It turns out that the endpoints A = or A = 1 are in general extensively degenerate 
eigenvalues of eq. (3.10). These eigenvalues do not correspond to resonances. Indeed, we 
know from section 2.1 that the conductance has no pole at A = 1, while it has a simple 
pole at A = if, and only if, the Q-phase is conducting. Let us set £q = 1 in this situation, 
£0 = else. The number of resonances then reads 

nR = ns -no -ni + Eq, (3.11) 



where no and ni denote the respective multiplicities of the endpoint eigenvalues A = 
and A = 1. More precisely, no is the number of zero- modes of the operator Aq, i.e., the 
dimension of its kernel. Since Aq has one zero-mode per cluster C of the Q-phase which 
is disconnected from the electrodes, we have 

no^ns-Y.{s{C)-l + x{C)) . (3.12) 

CCQ 

In this formula, the sum nms over the clusters C of the Q-phase, s{C) denotes the number 
of sites of the cluster C, and the characteristic function x(C) is unity if the cluster C 
overlaps with either of the electrodes, and zero else. The multiplicity rii of the eigenvalue 
A = 1 can be expressed similarly, in terms of the clusters of the P-phase. 

In the thermodynamic limit, one can derive from eq. (3.12) expressions for the fractions 
of eigenvalues which are condensed at A = and A = 1. Indeed, the terms x(C') are 
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negligible, while the other terms can be expressed as functions of geometrical characteristics 
of the bond percolation problem [22] , namely 

l-2q + P{q) + nc{q), 

(3.13) 

l-2p + P{p) + nc{p), 

where ndp) is the mean number of finite clusters per site, while P{p), the percolation 
probability, is the probability for any given bond to belong to the infinite cluster. The 
latter quantity is non-vanishing only for p > pc- 

The density of resonances pnip), defined in eq. (2.32), then reads 

Pr{p) = 1 - po{p) - pi(p) = 1 - P{p) - P{q) - n,{p) - n,{q)- (3-14) 

This quantity is symmetric in the exchange p g, in agreement with eq. (2.37). For a 
small concentration p, the contributions to eq. (3.12) of all the P-clusters and Q-clusters 
consisting of up to four bonds can be enumerated by hand. We thus obtain 

Po{p) =/ + ■■■, Piip) = 1-2P + / + ---, p^(p) =2p-2/ + --- (3.15) 

At the percolation threshold, Pr{p) takes its maximal value, which can be determined as 
follows. The percolation probability P{pc) vanishes, while ricipc) is known exactly [23], 
hence 

Po(Pc) = Pi(Pc) = ne(pc) = = 0.098076, pfl(pe) = 3(2 - Vs) = 0.803848. 

(3.16) 

3.3 Residues of the conductance 

At any resonance corresponding to a non-trivial eigenvalue (A^ 7^ and 1) of eq. (3.10), 
the map of the electric potentials on the network is given by the components (Xa)x of the 
associated right eigenvector X^. Since eq. (3.10) is symmetric, the are simultaneously 
its left and its right eigenvectors: 

X^(Ag-A„A) = 0, (Ag - A„A)X„ = 0, (3.17) 

with the row vector X* being the transposed of the column vector X. The ns eigenvectors 
{Xa} form a basis. They are orthogonal to each other with respect to the metric (—A), 
namely X* (— A)Xb = for a 7^ 6. We normalise them as 

X* (-A)X, = 5a,h. (3.18) 

The squared norm of the eigenvector X^ reads 

X^(-A)X„ = ^ = 1, (3.19) 
(x,y) 



Poip) 
Pi{p) 



lim 



no 



M,N^oo MN 

hm — 

M,N^oo MN 
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where 

E^,y = - {Xa)y (3.20) 

is the local electric field on the bond (x, y). 

The solution V to the inhomogeneous Kirchhoff equations (3.1), (3.6) can be obtained 
in terms of the eigenvectors Xa. Indeed, let us expand V on the basis of the X^: 



V = J]c,(A)X«. (3.21) 



a=l 



By inserting this expansion into eq. (3.6), and multiplying to the left by X^, we readily 
obtain the amplitude Cb{X) in the form 

(AA(^)-B(«))-X6 , ^ 



Finally, by inserting this result into eqs. (3.9), (3.7), we obtain the following formula for 
the residues eta > of the conductance at its non-trivial poles {Xa ^ and 1) 



M 

NXn 



a. = -^((A«aW - B(^)) . X,) ((A«A(«) - B(«)) • X,). (3.23) 



As recalled in the beginning of section 3.1, the conductance also exhibits a simple pole 
at A = 0, if the Q-phase is conducting. The corresponding residue ai, which yields the 
conductance of the Q-phase by means of eq. (2.3), is not directly given by eq. (3.23). It can, 
however, be determined from the other ones {a2, • • • , ct^^), by using the sum rule (2.5). 

4 NUMERICAL RESULTS 

We have shown in section 3 that the full rational A-dependence of the conductance of 
a finite binary network can be expressed in terms of the generalised eigenvalues and eigen- 
vectors of eq. (3.10). We have implemented this algorithm numerically, using the IMSL 
routine EIGZS, in order to obtain numerical data concerning several quantities of interest, 
which will be discuss successively throughout this section. The CPU time for solving the 
spectral problem for each sample grows rapidly with the system size, proportionally to n^, 
i.e., to for a square sample of size N x N. 

The optimal use of our algorithm therefore consists in obtaining good statistics on 
samples of moderate sizes. We have commonly used sample sizes such as = 16 or 
N = 20, and a statistical ensemble of several times lO'* samples, having some 10^ random 
bond conductances in total, a good enough statistics to obtain very accurate data. No 
observable systematic finite-size effects have been found, even in the critical regions {p — > pc 
and A^OorA— >^1). This observation is in agreement with the argument, given at the 
end of section 4.1, showing that the critical regions are indeed very small in the variable 
A. 
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4.1 Density of resonances and spectral function 

We have evaluated the spectral density of resonances p(p, x) and the spectral function 
H{p,x), by means of their respective definitions (2.31) and (2.26). The spectral function, 
which plays a central role in the Bergman-Milton formalism [4, 5], had only been the 
subject of a limited amount of work. In ref [7] it has been extracted from the imaginary 
part of the conductance, measured at a small but finite distance e from the negative real 
h-QxSs. The present method yields a direct measurement of the spectral function, avoiding 
especially any contamination from the delta-function at a; = [see eq. (2.27)], which can 
be discarded in an exact way. 

Figures 3 and 4 respectively show histogram plots of p(p, x) and xH{j>, x), for values of 
the concentration p ranging from 0.05 to 0.5. Indeed, the symmetry relations (2.33), (2.37) 
allow to restrict our attention to p < = 1/2. Each plot contain the accumulated data of 
an ensemble of over 20,000 configurations of a network of size = 16 x 16, corresponding 
to 10^ random bonds in total. The density of resonances p{p,x) exhibits the expected 
symmetry (2.37) under the transformation x <->^ 1 — x within a good accuracy, for all 
values of p. This demonstrates that we have used a large enough statistical ensemble of 
random networks. The product xH{p, x) does not possess such a symmetry (except for 
p = Pc — 1/2), while the EMA prediction (2.53), shown as a semi-circle in Figure 4, is 
symmetric under x <-> 1—x. Both p(p, x) and xH(p, x) exhibit a rich structure, down to the 
scale of the resolution (each plot contains 100 bins). It will be demonstrated more clearly 
in section 4.2 that they are non-vanishing over the whole allowed spectrum (0 < a; < 1). 

For a small enough concentration, the most salient structures in p{p, x) and H{p, x) can 
be predicted from the analysis of the resonances of isolated finite clusters. For consistency 
with ref. [2], we shall consider the regime q = 1 —p 0. To leading order in q, the relevant 
configuration consists of one isolated Q-bond embedded in a host lattice consisting of 
P-bonds. This one-bond cluster, shown as configuration A in Figure 5, has one single 
resonance, at A = Aai = 1/2, yielding the observed leading peak in p{p,x) and xH{p,x). 
This one-bond cluster has two possible orientations, but only the horizontal case yields a 
non- vanishing residue [2], to leading order as g — > 0. The spectral function, the density of 
resonances, and the mean number of resonances therefore behave as 

H{p,x)^p6{x- 1/2), p{p,x) ^2p6{x- 1/2) (p^O), (4.1) 

in agreement with the sum rules (2.49) and with eq. (3.15), to leading order in p. To second 
order in p, p{p, x) and H{p, x) consist of a countable infinity of discrete components (delta- 
functions), corresponding to the resonances of configurations consisting of two bonds, in 
arbitrary relative position and orientation [2], and so on. The most salient subleading 
peaks have been marked in Figures 3 and 4 (dashed verticals) by some of the resonances 
of the two-bond and three-bond configurations shown in Figure 5. The configurations A, 
B, and E are self-dual, while (C,D) and (F,G) form dual pairs. The resonances of these 
clusters, determined exactly along the lines of ref [2], are given in Table 1. 

Right at the percolation threshold {p = Pc = 1/2), the data for the spectral function 
(Figure 4e) agree with the exact analytical result (2.36), which coincides with the EMA 
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prediction. The accuracy of this agreement provides another check of the quahty of the 
numerical simulations. The density of resonances remains a non-trivial function p{pc, x) 
at the percolation threshold (Figure 3e). Integrating the data of this Figure leads to the 
estimate Pr{Pc) ~ 0.80, again in good agreement with the exact result (3.16). It is worth 
noticing that the data shown in Figures 3e and 4e are practically not affected by the critical 
singularities at the endpoints of the spectrum. Indeed, the size of the critical region, given 
by eq. (2.46), can be estimated to be a few times 10~^, i.e., smaller than the width of the 
first or last bin. 

As the concentration varies from p = Q io p = Pc = 1/2, the profiles of the density of 
resonances and of the spectral function deform in a progressive way. They get smoother 
and smoother, with their maxima moving in a continuous way. The dashed verticals are 
shown as guides for the eye on all plots of Figures 3 and 4, although they only label the 
most salient structures for a small enough concentration. The spectral function also gets 
progressively in better and better agreement with the EMA prediction. 

4.2 Lifshitz tails 

We have also investigated the behaviour of the density of resonances near the endpoints 
a; = and a; = 1, in order to check the prediction (2.63) of the Lifshitz tail. We have chosen 
to investigate the density of resonances p{p,x), rather than the spectral function, because 
p{p, x) can be expected on general grounds to exhibit a clearer signal, in analogy with the 
one-dimensional situation [19]. Furthermore, the statistics can be doubled by using the 
symmetry relations (2.37), and each sample requires less CPU time, since the calculation 
of the eigenvectors Xa is not required for p{p, x). 

Figure 6 shows a logarithmic plot of the integrated spectral density of resonances 

Pint{p,x)^ p{p,y)dy (4.2) 

^0 

against for a concentration p = 0.1. The data correspond to over 103,000 samples of 

size AT = 16 X 16, i.e., to 5 x 10^ random bonds. The range of the plotted data corresponds 
to a; < Amin and Amax ^ X < 1, with Amin and Amax being the endpoints (2.52) of the 
EMA prediction for the spectral function. Indeed, the Lifshitz behaviour is expected to 
manifest itself mostly out of the "bulk" of the spectrum, the latter being conveniently 
defined as the support of the EMA formula [20]. A linear behaviour is clearly observed, 
confirming the analytical form (2.63) of the Lifshitz tail. A further qualitative confirmation 
of Hesselbo's argument on the Lifshitz behaviour is as follows. The data of Figure 6 exhibit 
oscillations around to the fitted straight line, and the top of each of the most prominent of 
these oscillations corresponds, with a good accuracy, to the lowest resonance of the hairpin 
structures shown in Figure 2a, embedded in the square lattice, with n = 1 to 4 cells. The 
case n = 1 corresponds to configuration G of Figure 5, with its lowest resonance Agi- 

From a quantitative viewpoint, the slope of the fitted straight line in Figure 6 yields 
C{p = 0.1) 3.19. The amplitude C{p) has been similarly measured for p = 0.05, 0.15, 
and 0.2. The results are plotted in Figure 7 against |lnp|. The data are nicely fit to 
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the straight hne C{p) ^ 1.98|ln7?| — 1.39, to be compared with the analytical predic- 
tion (2.64). The slope 1.98 is some 10% smaller than the analytical value 11/^/2 = 2.2214. 
The intercept yields the estimate S ~ 1.40 for the configurational entropy per cell, a 
significantly larger value than the entropy of self-avoiding walks on the square lattice, 
'S'sAW = In //SAW = 0.970 [24]. These observations suggest that other types of linear ex- 
tended structures, besides the worm-like ones identified in the framework of Hesselbo's 
argument, may contribute to the Lifshitz behaviour of the conductivity. 

4.3 Distribution of spacings between resonances 

The distribution of spacings between successive energy levels has been extensively 
investigated in a variety of quantum systems, ranging from nuclei to billiards [25]. Generic 
spectra belong to three universality classes of level spacing distributions, according to their 
symmetry properties, in correspondence with the classical ensembles of real symmetric, 
Hermitian, and symplectic random matrices [26] , respectively called GOE, GUE, and GSE. 

We have investigated the distribution of spacings between successive resonances in the 
range 

p = p^ = l/2, l/4<x<3/4, (4.3) 

where the spectral density of resonances is very fiat, i.e., very close to being a constant, 
p = 0.658 (see Figure 3e). The range (4.3) can be considered as fully generic, although 
the concentration assumes its critical value Pc, since the critical singularities only influence 
very small regions around the endpoints a; = or a; = 1, as explained in section 4.1. 

For a finite network of size NxN, the mean spacing between two successive resonances 
is approximately Xa+i — Aa ~ l/{pN'^). We thus define the reduced spacings 

Sa = pN\K+l-Xa). (4.4) 

Figure 8 shows a histogram plot of the distribution P{s) of the spacings s^j, obtained from 
an ensemble of networks of size 20 x 20 having 10^ random bonds in the range (4.3) in 
total. This distribution obeys by construction P(s)ds = sP(s)ds = 1. 

Since the generalised eigenvalue problem (3.10) involves two real symmetric matrices, 
A and Ag, the natural universality class to which the data for P{s) should be compared 
is the Gaussian orthogonal ensemble (GOE). The law -Pgoe('S) [26] is shown as a full 
line in Figure 8. The data share characteristic qualitative features of the GOE spacing 
distribution, including a linear repulsion at short spacings, i.e., P{s) ~ s for s <^ 1, and 
a fast fall-off at large spacings. There is, however, a small but significant quantitative 
difference between the data in the range (4.3) and the GOE prediction. 

4.4 Number of visible resonances 

We now turn to the number of visible resonances of the RL — C model on a finite 
network. A resonance is said to be visible if it corresponds to a true maximum in the 
frequency dependence of the real part of the admittance, as given by eq. (2.20a). 

For a finite network of size N x N, the typical spacing between resonances scales as 
l/A?"^, as already mentioned in section 4.3. Since two resonances of comparable strengths 
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are resolved, i.e., separately visible, if their spacing is larger than the width 1/Q of each of 
them, we are led to propose the following finite-size scaling law for the fraction of visible 
resonances 

(^fl)visib _ p / Q 1 ^4 ^^ 



riR 

The scaling function F{X) is expected to grow monotonically from F(0) = 0, since a 
vanishing fraction of the resonances is visible if Q does not scale as A?"^, to F{oo) = 1, 
since all the resonances of a finite sample are eventually visible for a large enough quality 
factor. 

Figure 9 shows a plot of the fraction of visible resonances over the whole spectrum 
(0 < A < 1), for p — = 1/2. The observed collapse of the data for the sizes 12 x 12 
and 20 x 20 clearly demonstrates the validity of the scaling law (4.5). The full line shows 

a numerical fit of both series of data to the common analytical form l/(l — F{X)^ = 
1 + XP2{X), with P2{X) being a quadratic polynomial. The quality of the fitted curve, 
meant as a guide to the eye, suggests a linear behaviour of the scaling function at small 
X, as well as a power-law convergence of the form 

1 - F{X) ~ X-^/^ {X > 1). (4.6) 

The fraction of visible resonances depends a priori on how uniformly the resonances 
are distributed, and on the dispersion in the corresponding residues ««■ Both features 
can depend quantitatively on the concentration p. The scaling function F(X) is therefore 
expected not to be universal, but rather to weakly depend on the concentration p and on 
the range of values of A under consideration. Its main qualitative features, such as the 
power law (4.6), are however expected to be universal. The same remarks apply to the 
distribution of spacings between resonances, investigated in section 4.3. 

4.5 Distribution of local electric fields 

The algorithm presented in section 3 also gives access to the spatial structure of the 
resonances. Indeed, the eigenvector directly provides a map of the electric potentials 
at the resonance corresponding to the eigenvalue A^. For each resonance, we define the 
local electric field on the bond (x, y) as 

E^,y = V^-Vy^ C((X„)^ - (Xa)y) . (4.7) 

The electric fields are defined up to an overall multiplicative constant c. We choose this 
constant to be c = ^/nB■, so that the normalisation (3.19) of the eigenvectors implies 

E ^x,y = nB. (4.8) 
(x,y) 

One could think of many ways of analysing the spatial structure of the electric field at 
resonance. We have chosen to investigate the distribution of the local field living on any 
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given bond. More precisely, we have evaluated the successive moments of this distribution 
on square samples of size N x N, namely 

S,{N) = {\E\'^) = (-^J2\^->y\')- (4.9) 

The normalisation (4.8) ensures the identities So{N) = 5*2 (A^) = 1. 

Figure 10 shows a log-log plot of the first non-trivial moments, of index /c = 1, 3, 4, 5, 
and 6, against the linear size N of the sample, from AT = 6 to AT = 24. Data are obtained 
in the range (4.3), with around 10^ random bonds for each sample size. Power laws of the 
form 

Sk{N)^N'"' (4.10) 

are clearly apparent. This scaling behaviour, with a non-trivial dependence of the expo- 
nent Xk on the index k, is a signature of multifractality [27, 28]. Along the lines of the 
multifractal formalism, we introduce the generalised (Renyi) dimensions dk, such that 

X. = (1^3^. (4.11) 

The dk are expected to decrease from do = 2, the dimension of the network, to doo = 0. 

In physical terms, multifractality implies that the patterns of resonant electric fields 
exhibit strong local fluctuations, rather similar to those observed in wavefunctions of the 
Anderson model, in the marginal two-dimensional situation [29]. In particular the reso- 
nances are neither localised nor extended, in the conventional sense of these expressions. 
Indeed, extended patterns of electric fields would correspond to dk — 2 for all k > 0, while 
localised ones would have dk = ior k > 1. We recall that a similar phenomenon of 
multifractality has been reported for the DC problem of the conductor-insulator mixture, 
right at the percolation threshold [30, 31]. 

From a quantitative viewpoint, a more refined fit of the data to the power laws (4.10), 
including a relative correction in 1/A^, yields more accurate estimates for the exponents 
Xk and the associated dimensions dk, listed in Table 2. Figure 11 shows a plot of the dk 
against the index k. An approximate linear decay of the form 

dk^2-(3k, (4.12) 

with (3 ~ 0.194, is observed over a fairly broad range (0 < A; < 4). A similar linear 
behaviour has been predicted analytically for the two-dimensional Anderson model in 
the weak-disorder regime [29]. The linear law (4.12) corresponds to a log-normal (LN) 
distribution of the local fields. Indeed, let us set 

e = \ii\E\. (4.13) 

The scaling law (4.10) then reads (exp(A;£)) ~ exp{Bk{k - 2)/2), with B = j3lnN. Ne- 
glecting the A;-dependence of the prefactor, this last expression corresponds to a log-normal 
distribution for i.e., to a Gaussian law for the logarithmic variable £, namely 

nLN(^) = {2nB)-'/'exp (-^^^^) • (4-14) 
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The actual probability density n(£) is shown in Figure 12, for networks of size 16 x 16 
in the range (4.3). This very asymmetric distribution looks quite different from a log- 
normal law. In particular, it falls off as n(£) ~ exp(— as £ ^ — oo, yielding a finite value 
at \E\ — of the probability density P{\E\) — e^Il{£). Most of the dependence of the 
distribution on the sample size N takes place to the right of the plot, for large values of £, 
close to the upper bound £ma.x = (lnnB)/2, where a very fast decay is observed. 

5 DISCUSSION 

We have investigated the AC conductivity of binary random networks of complex 
impedances, with emphasis on its analytic structure in the complex variable h or A, and 
on the corresponding resonant behaviour. The present analysis is an extension to the 
general binary case of a previous work [2], devoted to the resonant response of a finite 
cluster, or set of clusters, embedded in an infinite homogeneous host lattice. Along the 
lines of refs. [1, 2], the poles of the conductance are interpreted in terms of the resonances 
which show up in the AC conductivity of the RL — C model, and of the relaxation times 
in the transient response of the R — C model. We have proposed an efficient algorithm, 
which allows a determination of the rational A-dependence of the conductance y of a finite 
sample, in terms of its poles A^ and of the associated residues Ua- A very similar formalism 
had been proposed long ago by Straley [6] . 

We want to underline again that the main advantage of the present approach is to give 
at once the analytic structure of the conductance in h or A, for any finite sample. As far as 
a numerical investigation of the resonant response is concerned, this approach is therefore 
more suitable than the usual numerical methods, which can only yield the conductivity of 
the binary model for a fixed value of the ratio /i, such as the transfer-matrix method [12], 
or the iterative algorithm based on the y — A transformation [32] , used e.g. in ref. [7] . It is 
also worth noticing that our approach yields the full analytic structure of the conductance, 
including the static conductance of the Q-phase of any finite sample. Indeed the latter 
quantity is given by eq. (2.3), with the corresponding residue cti being given by the sum 
rule (2.5) in terms of all the other residues, corresponding to genuine resonances. 

An extensive use of this algorithm, in the case of the square lattice, has allowed us to 
investigate in detail many aspects of the resonant dielectric response of the binary model. 
In general we have used 10^ random bonds or more per measurement, a good enough 
statistics to obtain very accurate data. We have investigated the density of resonances 
p{p,x) and the spectral function H{p,x). This approach yields a better evaluation of the 
spectral function than the more direct method, consisting in measuring the imaginary part 
of the conductivity. The most salient structures have been labeled, at least for a small 
enough concentration p, by resonances of configurations made of one to three bonds, which 
can be determined exactly, along the lines of ref. [2]. The data for the spectral function 
have also been compared with the EMA prediction. 

The Lifshitz behaviour of the density of resonances near the endpoints has been suc- 
cessfully characterised. A good qualitative agreement is found with the argument of Hes- 
selbo [21], according to whom the analogues of the Lifshitz sphere are linear extended 
objects, such as hairpins. From a quantitative viewpoint, our data suggest that other 
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classes of extended structures may also contribute to the Lifshitz tails, even for a small 
concentration p. The present situation is a lucky one, since numerical investigations of 
Lifshitz tails in more than one dimension are known to be a very difficult task in general, 
especially in the case of binary disorder [20]. The distribution of the spacings between 
neighbouring resonances has been found in qualitative agreement with the universal dis- 
tribution of the GOE universality class of random matrices, although small but definite 
differences show up at a quantitative level. The number of visible resonances of the RL — C 
model on a finite sample of size x A^, as a function of the quality factor Q, has been 
shown to obey a finite-size scaling law (4.5), involving a scaling function F[X) of the vari- 
able X = Q/N^. Quantities such as the scaling function F{X), or the spacing distribution 
P{s), are expected not to be universal: quantitative features of these functions should 
rather weakly depend on the concentration p and on the range of resonances considered, 
the range (4.3) being meant as a generic example. 

More generally, for the binary model on a (i-dimensional lattice, the appropriate finite- 
size scaling variable describing dielectric resonances reads X — Q/N'^. Indeed, the number 
of resonances on a sample of linear size N scales as its volume N^. This observation 
yields in particular a prediction for the divergence of the current correlation length ^ in 
the weak-dissipation regime (Q » 1). By setting X ~ 1 for A?" ~ ^, we get 

^~g^ with J^=^- (5.1) 

We thus recover a simple result due to Hesselbo [21], which has been corroborated by 
numerical simulations, yielding z7 = 0.4 ±0.1 in two dimensions [8]. 

Finally, we have investigated the distribution of the resonant electric field living on any 
given bond. The moments of this distribution obey power laws with non-trivial exponents 
Xk, a characteristic feature of multifractality. The associated generalised dimensions dk 
are found to behave similarly to those observed in the Anderson model of localisation, in 
the marginal two-dimensional case [29]. Multifractality thus appears to be a quite generic 
feature of the resonant response of binary networks. In particular, this phenomenon is 
unrelated to the percolation transition. In analogy with the fraction of visible resonances 
or the spacing distribution P{s), the exponents Xk and the dimensions dk are expected 
not to be fully universal, but to exhibit a weak dependence on the concentration p, and 
possibly on the range of resonances considered. 

This multifractal picture provides a quantitative characterisation of local features 
of the fluctuations in electric fields at resonance observed previously [8]. These giant 
fiuctuations have been argued to be responsible for surface-enhanced Raman scattering, 
this phenomenon being especially pronounced in strongly disordered semicontinuous films. 
There is a regime where the enhancement factor is predicted to be proportional to (E^) [9], 
whence the relevance of the dimension d^, in our notation. In a realistic system, with a 
small but finite dissipation rate 1/Q, these fiuctuations are expected to be critical, i.e., 
to exhibit strong spatial correlations, on scales smaller than the current correlation length 
^, estimated in eq. (5.1). Since the algorithm used in this work gives direct access to the 
full map of electric fields at resonance, it could be used to investigate other aspects of 
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dielectric resonances, including their spatial correlations, on which some information is 
already available [8]. 
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FIGURE CAPTIONS 



Figure 1: Schema of the binary network under consideration, with M = N = 8. The 
conductance Y is measured between the plane electrodes (bus bars). P-bonds with con- 
ductance (To are shown as solid lines, Q-bonds with conductance ai are shown as dotted 
lines. 

Figure 2: (a) hairpin configuration and (b) worm-like configuration, with n = 8 cells. 
Same conventions for the bond conductances as in Figure 1. The transversal conductance 
Yn is measured between the point electrodes shown as large dots. 

Figure 3: Histogram plot of the spectral density of resonances p{p,x), for (a) p — 0.05, 
(b) p = 0.1, (c) p = 0.2, (d) p = 0.3, (e) p = 0.5. The dashed verticals show some of the 
resonances, listed in Table 1, of the configurations made of one to three bonds, shown in 
Figure 5. 

Figure 4: Histogram plot of xH{p,x), with H{p,x) being the spectral function of the 
conductivity, for (a) p = 0.05, (b) p — 0.1, (c) p — 0.2, (d) p = 0.3, (e) p = 0.5. Semi- 
circles show the EMA prediction (2.53). Dashed verticals as in Figure 3. 

Figure 5: Configurations made of one to three bonds, whose resonances are listed in Table 
1, and used to label the most salient peaks in Figures 3 and 4, for p small enough. 

Figure 6: Logarithmic plot of the integrated density of resonances pintip,x), against 
for a concentration p = 0.1. The slope of the least-squares fit yields the amplitude 
C{p — 0.1) ~ 3.19. Numbers label the lowest resonances of hairpin structures, as explained 
in the text. 

Figure 7: Plot of the measured amplitude C{p) of the Lifshitz tail, against |lnp|. The 
fitted straight line is discussed in the text. 

Figure 8: Histogram plot of the distribution P{s) of normalised spacings between neigh- 
bouring resonances, measured in the range (4.3). The full line shows the GOE prediction. 

Figure 9: Plot of the finite-size scaling function F{X) of the fraction of visible resonances, 
against the scaling variable X = Q/N'^. Data correspond to samples of size N x N, with 
iV = 12 and = 20, at p = 1/2. The full line, showing the fit described in the text, is 
meant as a guide for the eye. 

Figure 10: Log-log plot of the moments Sk{N) of the local electric field at resonance, 
measured in the range (4.3), against the linear sample size N {6 < N < 24). 

Figure 11: Plot of the generalised dimensions dk of the distribution of electric fields at 
resonance, measured in the range (4.3), against index k. The straight line shows the linear 
behaviour (4.12), with P = 0.194. 

Figure 12: Logarithmic histogram plot of the probability density n(£) of £ = In the 
logarithm of the electric field at resonance, for samples of size 16 x 16, measured in the 
range (4.3). The dashed vertical marks the upper bound ^max = (lnnB)/2, with ub = 481. 
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TABLE CAPTIONS 



Table 1: Location of the resonances of some configurations consisting of one to three 
bonds, shown in Figure 5, and used to label the most salient peaks in Figures 3 and 4, for 
p small enough. 

Table 2: Exponents Xk and associated generalised dimensions dk characterising the mul- 
tifractal distribution of local electric fields at resonance, measured in the range (4.3). 
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Table 1 



configuration 


resonances 


A 


Aai 


= 1/2 = 0.500000 


B 


Abi 

Ab2 


= I/tt = 0.318310 
= 1 - I/tt = 0.681690 


C 


Aci 

Ac2 


= 1 - 2/7r = 0.363380 
= 2/7r = 0.636620 


D 


Adi 

Ad2 


= Aci 

= Ac2 


E 


Aei 

Ae2 
Ae3 


= 1/2 - \/2(l/2 - I/tt) = 0.243051 
— Aai 

= 1/2 + v^(l/2 - I/tt) = 0.756949 


F 


Api 

Af2 
Af3 


= 1/4 + I/tt - 'u;/4 = 0.302436 
= Aci 

= 1/4 + I/tt + w;/4 = 0.834184 


G 


Agi 

Ag2 
Ag3 


= 3/4 - I/tt - w/4: = 0.165816 

— Ac2 

= 3/4 - I/tt + w/4 = 0.697564 




w = 


v/9-40/7r + 48/7r2 
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Table 2 



7 

k 


exponent Xk 


dimension dk 








2 


-1 

1 


—0.097 ± 0.016 


1.806 ± 0.032 


2 







3 


r\ c\f\^ 1 r\ r\'\ o 

0.295 ± 0.018 


1.410 ± 0.036 


4 


0.7d9 ± O.OZD 


l.zol ± O.OZD 


5 


1.390 ±0.026 


1.074 ±0.018 


6 


2.130 ±0.024 


0.935 ±0.012 
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